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ABSTRACT 

o\ _ 

The first results from a Tully-Fisher (TF) survey of cluster galaxies are pre- 
sented. The galaxies are drawn from fifteen Abell clusters that lie in the redshift 
| 7 | ■ range ~ 9000-12,000 kms -1 and are distributed uniformly around the celestial 

sky. The data set consists of i?-band CCD photometry and long-slit Ha spec- 
troscopy. The rotation curves (RCs) are characterized by two parameters, a 
CN ■ turnover radius r t and an asymptotic velocity v a , while the surface brightness 

profiles are characterized in terms of an effective exponential surface brightness 
I e and scale length r e . The TF scatter is minimized when the rotation velocity is 
. measured at r = (2.0±0.2)r e ; a significantly larger scatter results when the rota- 

tion velocity is evaluated at ^,3 or £ 1.5 scale lengths. This effect demonstrates 
CT\ ■ that RCs do not have a universal form, specified only by luminosity, as suggested 

by Persic, Salucci, & Stel. The scatter minimum at r = 2r e is interpreted in 
terms of a thin stellar disk plus spherical dark halo mass model. Variations in 
halo mass and size at fixed disk mass and size can produce extra TF scatter at 

■4-^ i 

arbitrary radii, but no additional scatter at the special radius r = 2r e provided 



Ph oc M^ ' 7 . In contrast to previous studies, a modest but statistically signifi- 
cant surface-brightness dependence of the TF relation is found, v oc L°' 28 I® AA . 
This indicates a stronger parallel between the TF relation and the correspond- 
ing Fundamental Plane relations of elliptical galaxies than has previously been 
recognized. The scatter of the optimized TF relation decreases with increasing 
luminosity and surface brightness, from ~0.75 mag for low-luminosity, low-SB 
objects to ^0.35 mag for high-luminosity, high-SB objects. This effect is well 
described by a model in which the intrinsic TF scatter is ~ 0.30 mag, and most 
of the remaining scatter is caused by rotation velocity measurement errors of 
~ 15 kms -1 independent of rotation amplitude. Future papers in this series 
will consider the implications of this cluster sample for deviations from uniform 
Hubble flow on £ 100/iT 1 Mpc scales. 
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1. Introduction 

In the late 1980s and early 1990s, several groups reported the existence of large-scale 
bulk flows in the local universe (Dressier et al. 1987; Willick 1990; Mathewson, Ford, & 
Buchhorn 1992; Han & Mould 1992; Courteau 1992; Courteau et al. 1993; Lauer & Postman 
1994; recent reviews of the subject include Dekel 1994, Postman 1995, Strauss & Willick 1995, 
and Strauss 1997). These studies suggested that high-amplitude (v p ^,500 kms" 1 ) peculiar 
velocities might be coherent across large volumes, tens to perhaps hundreds of megaparsecs 
in diameter. The actual coherence scale of the flows has been poorly constrained, however. 
The results of Lauer & Postman, in particular, suggested that this scale could be as large as 
300/T 1 Mpc. Other, more recent, studies (e.g. Riess et al. 1997; Giovanelli et al. 1998a,b) 
have indicated a much more modest coherence scale for bulk flows, ^ 50/i -1 Mpc. 

The question of the coherence scale for high-amplitude bulk flows is an important one 
for cosmology, because this scale is sensitive to the underlying power spectrum of density 
fluctuations. In most models of large-scale structure formation, one does not expect Lauer- 
Postman like flows (see, e.g., Strauss et al. 1995). To address the question of the bulk flow 
coherence scale, the author undertook a Tully-Fisher (TF) and Fundamental Plane (FP) 
survey of cluster galaxies beginning in 1992. The observations for this survey (§3) were 
carried out from the Las Campanas (LCO) and Palomar Observatories, and the cluster 
sample was restricted to a narrow range of redshifts near cz = 10,000 kms -1 . For this 
reason, the survey is referred to as the LCO/Palomar 10,000 km s _1 Cluster Survey (LP10K). 
Observations and analysis of the TF data have recently been completed, and the first results 
are reported here. A second paper (Willick 1998a, Paper II) will discuss constraints on bulk 
flows from the LP10K TF data set, and a third paper will present the TF data (Willick 
1998b, Paper III). The FP survey is ongoing, with results expected in 2-3 years. 

The focus of the present paper is on the form of the TF relation that emerges from 
the LP10K data set. Originally, TF studies used photoelectric apparent magnitudes and 
21 cm velocity widths. Thus, the amount of data available for each galaxy was limited. 
The LP10K TF survey is one of a number of recent TF surveys in which CCD surface 
photometry is combined with long-slit optical spectroscopy. The information per galaxy is 
thus considerably increased relative to the original TF surveys, and the implications of these 
additional data are only now being fully appreciated. A key goal of this paper is to explore 
how to make optimal use of the full range of photometric and spectroscopic data in applying 
the TF relation. 

The outline of this paper is as follows. § 2 provides a brief overview of the key issues 
involved in analyzing TF data, and discusses the significance of these issues for our under- 
standing of galaxy structure. §§ 3 and 4 describe the observations and data reduction for 
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the LP10K survey. § 5 presents the TF fitting procedure and results. § 6 presents an inter- 
pretation of the key features of the TF relation in terms of galaxy structure. § 7 summarizes 
the main points of the paper. 

2. The TF Relation 

The relationship between spiral galaxy luminosity and rotation velocity (the TF relation) 
was first noted by Tully & Fisher (1977). They and subsequent workers in the 1970s and 
1980s measured rotation velocity from the width of the 21 cm line (see, e.g., Aaronson et 
al. 1982, 1986; Bottinelli et al. 1983; Pierce & Tully 1988). During this period few TF 
studies made explicit use of the fact that, in reality, spirals do not have a single, well-defined 
rotation velocity, but rather exhibit a rotation curve (RC), v c (r), where r is distance from 
the galaxy center. However, the fact that many spirals exhibit "flat" RCs (e.g., Rubin et al. 
1982), in which v c (r) rises rapidly from zero at the center and then is constant to the largest 
observed radii, suggested that the 21 cm width did indeed reflect a well-defined quantity: 
"the" rotation velocity — i.e., the velocity of the flat part of the RC — of the galaxy. 

In the late 1980s and early 1990s, several new TF surveys substituted, in whole or in 
part, long-slit optical spectroscopy of the Ha emission line for 21 cm velocity width measure- 
ment. Courteau (1992, 1997) analyzed over 300 Northern sky Sb-Sc galaxies using i?-band 
photometry and long-slit spectra. Courteau (1997; see also Courteau & Rix 1999) showed 
that this data set produced an optimized TF relation when the velocity was measured from 
the RC at 2.2 exponential disk scale lengths. Mathewson, Ford, & Buchorn (1992) carried 
out a TF study of over 1300 Southern sky objects, of which about half had optically mea- 
sured RCs. This subset was reanalyzed by Kasen (1997) who found an optimized TF relation 
when the rotation velocity was measured at 1.8 disk scale lengths. Schlegel (1995) conducted 
a TF survey of IRAS galaxies and found an optimized TF relation for v c (r) evaluated at 
an isophotal radius. Simon (1998) analyzed a 95-galaxy subset of the SHELLFLOW sample 
(Courteau et al. 1998) and found a minimum TF scatter when the the rotation velocity was 
measured at 1.7 disk scale lengths. 

While differing in detail, the above studies all demonstrated that there is no unique, a 
priori TF rotation velocity for spirals with carefully measured RCs. In particular, they agreed 
that fitting the observed RCs and using the asymptotic rotation velocity, the amplitude of the 
flat portion of the RC, does not yield an optimized TF relation. Rather, the RC amplitude 
must be evaluated at a radius determined from the photometric properties of the galaxy. 
A corollary is that even galaxies whose RCs never "turn over" — are still rising even at the 
outermost radii observed — can be assigned rotation velocities that fit the TF relation. The 
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analysis presented here will confirm and extend these conclusions. These issues are important 
not only for applying the TF relation, but also for understanding spiral galaxy structure, 
because the radius at which velocity best correlates with luminosity is a function of the 
relative distribution of dark and luminous matter (Courteau & Rix 1999; McGaugh & de 
Blok 1998; Navarro 1998). 

A second element of recent TF surveys has been accurate CCD surface photometry. 
Thus, not only apparent magnitudes, but well-defined surface brightnesses, radii, and lumi- 
nosity shape parameters are available for all galaxies. In principle, this photometric infor- 
mation allows the incorporation of additional parameters into the TF relation. However, 
efforts to find such parameters have not yielded fruit in the past (Han 1991; Willick 1991). 
In particular, the relation has been found to be independent of disk surface brightness (e.g., 
Tully & Verheijen 1997). The picture that has emerged has been of the TF relation as a two- 
parameter (one-dimensional) relation between rotation velocity and luminosity only. This 
picture contrasts with that which has emerged for elliptical galaxies over the last decade 
(e.g., Djorgovski & Davis 1987; Bender, Burstein & Faber 1992; Jorgensen et al. 1996). El- 
lipticals exhibit a three-parameter relation between luminous radius, surface brightness, and 
velocity dispersion known as the Fundamental Plane (FP). This has led to the question of 
why ellipticals and spirals differ in this regard. An important conclusion of this paper will be 
that spirals in fact exhibit properties more similar to ellipticals than previously suspected; 
when properly analyzed the TF relation is a three-parameter one similar to the FP. 



3. Observations and Preliminary Data Reductions 

The observations were carried out during the period 1992-1995 at the Las Campanas 
Observatory (LCO) in Chile and at Palomar Observatory in California. Fifteen Abell Clus- 
ters were selected as target fields from which the galaxy samples were to be selected. They 
are shown in Table 1. The selection criteria were (1) published cluster redshifts in the range 
0.03-0.04 (9000 < cz < 12, 000 kms" 1 ), (2) an approximately isotropic, full-sky distribution, 
and (3) Galactic latitude |6| > 20° . The last criterion was imposed to minimize the effects of 
Galactic extinction. In practice, a large number of Abell clusters with 0.03 < z < 0.04 was 
identified, and 15 were selected from among them in such a way as to achieve the greatest 
degree of isotropy. 

The observational strategy was as follows. Kron-Cousins i?-band CCD imaging of wide 
(~ 1-2 square degree) fields centered on each cluster was carried out from the LCO 1 m and 
Palomar 1.5 m telescopes. Approximately 25 contiguous frames, each consisting of 2 or 3 
exposures totalling 15 minutes (LCO) or 10 minutes (Palomar), were taken in each cluster. 
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TABLE 1 

LCO/PALOMAR CLUSTER SAMPLE 



Abell # 


RA (1950) 


DEC (1950) 


z 


R a 


m\ 


Las Campanas Clusters c : 










2731 





07.7 


-57 


16 


.0312 





15.3 


3202 


3 


59.0 


-53 


18 


.0388 


1 


15.8 


0496 


4 


31.3 


-13 


21 


.0320 


1 


15.3 


3381 


6 


08.1 


-33 


35 


.0382 


1 


14.7 


2063 


15 


20.6 


8 


49 


.0337 


1 


15.1 


1139 


10 


55.5 


1 


46 


.0383 





15.0 


3578 


13 


54.7 


-24 


29 


.0400 


1 


15.1 


3733 


20 


59.0 


-28 


15 


.0386 


1 


15.6 


3869 


22 


18.2 


-55 


23 


.0396 





15.3 


2657 


23 


42.3 


8 


52 


.0400 


1 


14.9 


Palomar 


Clusters^: 












0260 


1 


49.0 


32 


55 


.0348 


1 


15.8 


0576 


7 


17.3 


55 


50 


.0381 


1 


14.4 


1228 


11 


18.8 


34 


36 


.0350 


1 


13.8 


2199 


16 


26.9 


39 


38 


.0300 


2 


13.9 


2247 


16 


52.0 


81 


39 


.0384 





15.3 



Table 1: Notes: (a) Abell Richness class, (b) Photographic magnitude of the 10th brightest 
cluster member, (c) Clusters observed primarily or exclusively from LCO. (d) Clusters 
observed primarily or exclusively from Palomar. 

The CCD frames were bias-corrected, flatfielded, registered and coadded using standard 
procedures within IRAFQ. Automated galaxy identification using the FOCAS package ( Jarvis 
& Tyson 1981) was then done, and all objects brighter than rrtR ~ 17 and minor-to-major axis 
ratios b/a ^ 0.2 were visually inspected. From these a subsample whose appearance indicated 
relatively late-type ( ^ Sb) morphology was selected for follow-up long-slit spectroscopy. The 
morphology cut was found to be critical for maintaining a high Ha detection rate. The precise 
magnitude limit varied from cluster to cluster, because an effort was made to obtain 15-20 



3 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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TF data points per cluster, which occasionally required observations of galaxies as faint as 
mji ~ 18. In addition to the cluster images, standard star fields from the compilation of 
Landolt (1992) were multiply observed on each LCO 1 m and Palomar 1.5 m night deemed 
photometric, providing zero points for the R band galaxy photometry. 

The long-slit spectroscopy was conducted from the LCO 2.5 m and Palomar 5 m tele- 
scopes at a resolution of ~ 1 A/pixel. The slit was oriented along the major axis of the 
galaxy as determined from the CCD images. The wavelength coverage was ~ 6200-7200 
A, allowing detection of Ha out to a redshift of z ~ 0.1. Exposure times ranged from 10- 
45 minutes (see below). The 2-dimensional spectrograms were flat-fielded, rectified, and 
wavelength-calibrated using standard IRAF procedures. Measurement of a rotation curve 
(RC) was possible only when extended emission lines were detected, about 50-60% of all 
exposures. Another ~ 10% of exposures yielded nuclear emission that sufficed for redshift 
determination but not RC measurement. In virtually all cases Ha was the strongest emission 
line seen and was the only one used to trace the RC. In a handful of cases a cosmic ray or 
chip flaw compromised Ha and made it necessary to use the N[II] or S [II] lines to trace the 
RC. 

Two unexpected factors complicated the observational program. First, the rate of de- 
tection of Ha was considerably lower than anticipated. To minimize waste of large-telescope 
time, a strategy was adopted in which short exposures (5-10 minutes) were taken, the result 
inspected for Ha, and a longer frame taken only if Ha was detected in the short frame and 
deemed likely to produce a usable RC. The relatively low Ha detection rate also meant that 
more candidates than expected had to be observed in order to obtain the hoped-for number 
of detections. Accordingly, the long exposures themselves, initially 45 minutes in length, 
were shortened to 20-30 minutes beginning in 1993. Moreover, the Ha detection rate varied 
significantly from cluster to cluster. The original plan called for ~ 15 TF objects per cluster. 
In practice, the number varies from as few as 8 objects to as many as 23 objects with TF 
data in the different clusters. These differences reflect the fact that in some clusters the Ha 
detection rate was over 90%, while in others it was ,$,30%. 

The second unanticipated complication was the prevalence of background objects in 
many of the clusters. The reason for drawing the galaxy sample from cluster fields was 
to ensure that a narrow redshift range around 10,000 kms" 1 was sampled. It turned out, 
however, that in some of these fields many galaxies did not lie at the published cluster 
redshift. In Abell 3733 (published redshift cz = 11,600 kms" 1 ), for example, the majority 
of the final TF sample objects have cz > 20, 000 kms -1 . Similar though less extreme results 
were found in a number of other clusters, especially those in the Southern celestial sky (the 
reason for this is unknown and may be worth further investigation). In the final TF sample 
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only about two-thirds of the objects have redshifts 7000 < cz < 15, 000 kms -1 , which might 
be called the extended target range. The remaining 1/3 almost all have 15,000 < cz < 
30,000 kms -1 , with only a handful of foreground (cz < 7000 kms -1 ) galaxies. Although 
this fact reduces the sensitivity of the survey to the bulk flow of a 10,000 kms -1 shell, it 
does not compromise the purposes of the present paper which are simply to explore the TF 
relation. 



4. Production of the Final Data Set 

4.1. Derivation of Photometric Parameters 

The images of galaxies with emission lines in the long-slit spectra were subjected to 
further analysis using the VISTA image-processing package (Stover 1988; see Willick 1991 or 
Courteau 1992 for a detailed discussion of spiral galaxy surface photometry using VISTA). 
Elliptical isophotes were fitted to the high to moderate S/N portion of the galaxy images. 
The fitted ellipses which best matched the galaxy disk were visually identified and were 
taken to define the global ellipticity (e = 1 — b/a) and position angle (PA) of the galaxy. 
The surface brightness profile was then extended to a final radius, rj, within ellipses of this 
fixed shape and orientation. Aperture magnitudes were also computed, via direct summation 
of pixel values, within these fixed ellipses out to rj. The value of rj was determined by the 
requirement that the galaxy surface brightness at this radius equal the uncertainty in the 
sky background determination. 

The surface brightness profile was extrapolated beyond r/ by fitting an exponential 
function, I[r) = J e _r//rd , to the outer portion of the measured profile. It was usually not 
appropriate to fit the entire profile because of manifestly nonexponential features such as a 
strong central bulge. The fit parameters Jo and r^, while not necessarily characteristic of the 
entire profile, were thus suitable for calculating the luminosity extrapolation, 



AL = 2ir(l-e)I r 2 d e 



-rf/r d 



1 + 2 

rd 



(1) 



which was added to the aperture magnitude within r j to obtain the "total" magnitude of the 
galaxy, later used in the TF analysis. Typical extrapolations via this method were 0.03-0.10 
mag. 

The fit parameters Io and r^, while suitable for the luminosity extrapolation, were not 
found to be the optimal measures of SB and radius for the TF analysis; as noted above, 
these parameters were often not characteristic of the bright inner parts of the galaxy. An 
alternative procedure for determining a characteristic surface brightness I e and scale r e was 
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therefore developed. The details are somewhat technical and are presented in Appendix A. 
Here we note only that (1) I e and r e are determined from the intensity moments of the galaxy, 
not from any kind of parametric fit, and thus are robust and objective, and (2) I e and r e are 
defined in such a way that they agree with the exponential disk parameters I® and if the 
entire surface brightness profile is well described by an exponential law. Thus, I e and r e are 
an "effective" exponential surface brightness and scale length, but they are not predicated on 
the galaxy's actually having an exponential profile. Their value as characteristic parameters 
is best justified a posteriori from their role in defining the TF relation (§5). 

A final photometric parameter that proved useful was the luminosity concentration index 
c, defined by 

c = 51og— (2) 

r 2 o 

where rx is the radius containing X% of the total light. 



4.2. Derivation of Spectroscopic Parameters 

The two-dimensional wavelength-calibrated spectrograms were interactively processed 
to extract the RC. Full details of the procedure are given by Simon (1998). Briefly, following 
sky-subtraction, the approximate position of the Ha emission was interactively marked out 
along the slit. Accurate central wavelengths were then found as a function of position, with 
spatial averaging done in regions of low S/N. The galaxy center was computed by finding 
the centroid of the continuum longward and shortward of Ha and averaging. A rotation 
curve was then obtained by applying the Doppler shift formula and subtracting the velocity 
corresponding to the position of the galaxy center, which was taken as the systemic redshift 
of the galaxy. 

The full RCs contain a significant amount of information, not all of which can be used 
in the TF analysis. It was thus decided to fit each RC with a simple function that would 
encompass its key features. A suitable choice is a two-parameter arctangent fit, 

v c {r) = 2 -^ tan- 1 (-) , (3) 
7r \r t J 

which was found to provide a reasonable fit to the RC in all casesQ. In Equation (|3|), v a 
represents the asymptotic or "flat" value of the RC, while r t is a "turnover" radius where 



In practice, a four-parameter fit was carried out, in which the galaxy center and the mean radial velocity 
were allowed to vary from their initial estimates obtained from the spectrogram processing. In general, the 
refinements of the central position and velocity were quite small. 
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the RC reaches one-half its asymptotic value, roughly separating the linear and flat part of 
the RC. It is important to note that v a is a fit parameter and not the actual asymptotic 
value of the RC, which, indeed, was not observed in many cases because the RC was still 
rising at the last observed radius. The arctangent fit is above all a convenient interpolation 
formula for determining the amplitude of the RC at any chosen radius, as discussed in the 
next section. 



4.3. Corrections to the Raw Data 

The procedures above yield the basic observational data for each galaxy: 

1. The kinematic data: heliocentric redshift cz & (kms -1 ) and RC parameters v a (kms -1 ) 
and r t (arcsec). 

2. The photometric data: total magnitude rriR (i?-band magnitudes) and effective surface 
brightness \x e (i?-band magnitudes per square arcsecond, the magnitude equivalent of 
I e ; cf. §4.1), characteristic radius r e (arcsec), concentration index c, and ellipticity 
e — 1 — b/a. 

Prior to fitting the TF relation several standard corrections were applied. The he- 
liocentric redshift cz Q was converted to a microwave background frame (CMB) redshift, 
henceforward denoted cz with no subscript. The galaxy inclination was computed from the 
ellipticity as described by Willick et al. (1997a). The asymptotic velocity v a was corrected 
for inclination and redshift broadening. The effective surface brightness ji e was corrected for 
(1 + z) A cosmological dimming. The total magnitude was corrected for Galactic extinction, 
internal extinction, and cosmological effects following Willick et al. (1997a), with two key 
differences. First, the coefficient of log(a/6) used in the internal extinction correction was 
1.20, somewhat higher than the value of 0.95 used by Willick et al. The higher value was 
chosen to eliminate any trends of the TF residuals with axial ratio, and may result from the 
fact that R band magnitudes have a shorter effective wavelength than the r band magni- 
tudes used by Willick et al. Second, the extinction maps of Schlegel et al. (1998), which are 
based on IRAS/DIRBE maps of dust emission, were used in preference to Burstein-Heiles 
extinctions. The corrected total magnitudes are denoted m with no subscript. 

An additional inclination correction was required for the surface brightness. The values 
of n e corrected only for cosmological dimming exhibited a significant correlation with axial 
ratio. This trend was removed by computing corrected effective surface brightnesses //^ = 
fi e — 0.85 log(a/6). The corrected values of fi e are used in the remainder of this paper. A 
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detailed justification of this and other corrections applied to the data will be presented in 
Paper III. 

5. Modelling the TF Relation 

5.1. Sample Definition 

The total LP10K TF sample, meaning all objects with both photometric and RC data, 
numbers 268 distinct galaxies. About fifty of these galaxies have multiple (2-4) photometric 
measurements, and about thirty have two RC measurements. (A detailed accounting will be 
provided in Paper III.) For the purposes of the TF analysis, we will consider as "objects" 
all photometry / spectroscopy data point pairs for a given galaxy. Thus, for example, if a 
galaxy has two RC measurements and two photometric ones, it contributes four objects to 
the analysis. The effect of this on the statistical significance of the results is discussed below. 

The multiplicity of photometry and spectroscopy is such that the 268 galaxies yield 386 
objects. To improve reliability we further pruned the sample as follows: (i) We required 
log(a/6) > 0.06 to exclude objects that are too face-on. (ii) We excluded objects with 
absolute magnitudes (based on a Hubble flow distance model, see below) fainter than —18.5 + 
5 log/i (galaxies fainter than this were consistently poor fits to the TF relation), (iii) Objects 
with purely linear RCs — for which the arctan fit constrained only the ratio v a /r t , not either 
parameter separately — were excluded, (iv) Objects whose RCs exhibited no perceptible 
width were excluded, (v) One galaxy with cz < 2000 kms^ 1 was excluded. These cuts 
eliminated 34 objects, leaving a full TF sample of 352 objects, comprising 245 distinct 
galaxies. 

The fact that a significant fraction of the sample has multiple measurements reduces 
the number of statistical degrees of freedom for the analysis. Normally, a sample of N data 
points and p free parameters has N—p degrees of freedom. This would still be the case here if 
the TF scatter arose entirely from observational errors, because then different measurements 
of an individual galaxy would be statistically independent. However, as will be shown below, 
a significant portion of the TF scatter is intrinsic. As a result, the TF errors of two or more 
data points corresponding to an individual galaxy are strongly correlated. 

A full discussion of this difficult issue will be deferred to Paper II. For the present we take 
a conservative approach and simply assume that the effective number of independent data 
points is simply the number of distinct galaxies in the sample, 245. The actual number is 
larger, because the observational scatter arises largely from the velocity width measurement, 
so we could arguably count the same galaxy twice if it has two measured RCs. Moreover, 
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the different photometric measurements of a given object, while they yield similar apparent 
magnitudes, produce somewhat different measures of inclination and scale length which yield 
different velocity widths even for a single RC. However, the conservative approach suffices 
here because the effects considered in this section have high statistical significance, and slight 
underestimates of their significance are thus unimportant. 

In what follows we use a likelihood analysis to study the properties of the TF relation. 
The likelihood statistic C we define in §5.4 has the property that, if the data points entering 
into it are independent, a decrease of one unit, over and above degrees of freedom added, 
corresponds to a one-sigma likelihood increase; an increase of four units over and above the 
added degrees of freedom corresponds to a two-sigma likelihood increase (cf. Willick et al. 
1997b for a more detailed discussion). To correct for the nonindependence of the data points, 
we scale the statistic by a factor N indep /N to t ~ 245/352 = 0.696. This scaling ensures that 
the above rules for assessing the significance of changes in C continue to hold, and we use 
these rules in the discussion to follow. 



5.2. Distance Assignments 

To fit the TF relation distances must be assigned to each galaxy from redshift-space 
information. The LP10K TF sample consists, in principle, of cluster galaxies, and the 
members of a given cluster are often assumed to lie at a common distance. However, for a 
variety of reasons this "cluster paradigm" (Willick et al. 1995) is a poor approximation here. 
As already noted, one-third of the sample consists of background galaxies with redshifts as 
high as z — 0.1. Applying the cluster paradigm would require a careful pruning of objects 
using distinct redshift limits for each cluster, which in the end would be rather arbitrary. 
Moreover, the true redshifts of the clusters are not necessarily well known. In addition, many 
sample spirals come from the outskirts of the clusters rather than from their presumably 
virialized cores. Finally, Willick et al. (1995) found that the redshift-distance relation for 
cluster spirals was often better approximated by Hubble flow than by the cluster paradigm 
in any case. 

For the above reasons, we assign distances (in Mpc) to all galaxies using a simple 
Hubble-flow model, 

d = H^cz , (4) 

where for definiteness we take H = 65 kms" 1 Mpc -1 . (Uncertainty in the Hubble constant 
translates directly into errors in the TF zero point derived from the analysis, but does not 
affect any other relevant parameters.) Modifications of this distance law for possible bulk 
flows, as well a careful division between cluster and non-cluster galaxies, will be considered 
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Fig. 1. — Two representative rotation curves from the LP10K sample. Both are from the cluster 
A2199. The left hand panel shows a galaxy with a classical flat RC, and with xt <C 1. The right 
hand panel shows a galaxy with a slowly rising RC, in which the asymptotic velocity is not reached 
in the observed region. Correspondingly, it has xt > 1. See text for further details. 

in Paper II. These considerations have no meaningful effect on the properties of the TF 
relation itself, however. 



5.3. Parameterization of the TF Relation 

Throughout this paper, we use the "inverse" form of the TF relation (Strauss & Willick 
1995, §6.1) in order to minimize selection bias effects. We define the "circular velocity 
parameter" 

77 = log (2t; TF ) - 2.5 , (5) 

where i>tf is the measure of rotation velocity that optimizes the TF relation (see below). 
The inverse TF relation is then written, in its simplest, two-parameter form, as 

r} = -e(M-D), (6) 

where M is absolute magnitude, and e and D are the inverse TF slope and zero point. 

We follow Courteau (1997) in supposing that Vtf is obtained by evaluating the fitted RC 
at a multiple of the characteristic radius of the galaxy. Here, however, we use the effective 
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exponential radius r e rather than the fitted exponential scale length (cf. §4.1). Thus, 

2f a _i / f s r e \ ,„s 

Vtf = tan , (7) 

vr V r t J 

where we treat f s as a free parameter to be determined by optimizing the TF relation. It 
will prove useful to define a dimensionless parameter 

x t = - , (8) 

a ratio of dynamical to luminous scale length. In terms of Xt we write the TF rotation 
velocity as 

VTF = ^ tan -i(k) . ( 9 ) 

7T \x t J 

The value of the x t parameter can best be appreciated from representative RC plots, as 
shown in Figure |IJ Because galaxy rotation is typically detected out to ~3-5 r e , the value 
of x t determines whether the curve exhibits a classical flat shape (x t < 1), or whether the 
RC is still rising at the outermost observed radii (x t ^ 1). Hence, x t is basically an RC shape 
parameter. We will make further use of this below. 



5.4. Method of Fit 

We adopt a maximum likelihood approach for determining the TF parameters. The 
method is based on the VELMOD likelihood approximation outlined by Willick & Strauss 
(1998; cf. Appendix 1 of that paper). The observed values of rj are assumed to be normally 
distributed about their predicted values, 

p( | x 1 / fa-[-e(m t -51og^-25-D]) 2 | 

P{rji\m u czi) = — == exp < — 2 \ . (10) 

V ^TCCTcS,i { Z<J eS,i ) 

In Eq. [10] er c fr is the effective TF scatter given by 

•^ 2 + G^f) 2 

where a v is the inverse TF scatter (intrinsic plus observational errors), and a v is the velocity 
dispersion relative to Hubble flow. We treat free parameter, but fix a v at 250 kms . 

The sample is too distant to reliably determine a v (cf. the discussion by Willick et al. 1997b). 
The likelihood for observing the entire data set is V = FL -P^l? 7 ^? c <2«)> where the product 
runs over all data points. We maximize this likelihood by minimizing C = —2 In V with 
respect to the various free parameters. 
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Fig. 2. — Inverse TF scatter (top panel) and the fit likelihood statistic C (bottom panel) for a two 
parameter (rotation velocity-luminosity) TF fit, plotted against log/ s , where f s is the parameter 
that determined the radius at which utf is evaluated (see equation ^ and the surrounding text for 
details). A strong minimum in TF scatter, and a corresponding maximum in fit likelihood, are seen 
at f s = 2.1. The dotted line in the lower panel shows indicates AC = 4.0 relative to the minimum 
of C, corresponding to the 2-a confidence limits on the value of f s . 



5.5. Establishing the value of f s 

The free parameters appearing explicitly in the model (equation [H]) at this point are D, 
e, and a v . In subsequent sections, we will treat f s , and other parameters not yet introduced, 
as free parameters. Before doing so, however, it is useful to maximize likelihood with respect 
to changes only in the zero point, slope, and scatter of the TF relation, for a range of 
fixed values of f s . This exercise will demonstrate the presence of a strong minimum in the 
TF scatter for f s « 2, a result that will prove robust to the introduction of addtional free 
parameters. 
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In Figure 0, the TF scatter a n and the corresponding value of C are plotted versus log f s 
for a sequence of fits done at fixed values of f s . Both a v and £ are minimized (likelihood is 
maximized) for f s = 2.0 ± 0.2 (one-sigma errors; the corresponding two-sigma, determined 
by the dotted line in the lower panel, are ~ ±0.35. Values of f s ^ 1.5 and f s ^ 3 are strongly 
ruled out. Note, in particular, that the TF scatter is much larger, and the likelihood of the 
fit vastly smaller, for very large (^10) f s , which corresponds to using the asymptotic flat 
portion of the RC as the TF rotation velocity. We conclude that there exists a characteristic 
radius, ~ 2 effective exponential scale lengths, at which rotation velocity correlates most 
strongly with luminosity. This conclusion is consistent with those reached by Courteau 
(1997), Kasen (1997), and Simon (1998). 

5.6. SB/Concentration Dependence of the TF Relation 

The maximum likelihood values of D, e, a v , and f s , for the two-parameter, rotation 
velocity-luminosity TF relation are given in the first row of Table 0. (Note that fitting the 
two-parameter TF relation involves four free parameters, because f s and a n also are varied.) 
In Figure [3] inverse TF residuals, ^(observed) — r\ (predicted), from this fit are plotted versus 
surface brightness (top panel) and concentration index (lower panel). A mild but significant 
trend is can be seen with respect to surface brightness. The diagonal line of slope 0.05 
indicates the approximate slope of this trend. A weaker trend is seen in the residual versus 
concentration index plot, though it appears to be significant only for the most concentrated 
galaxies (c^, 3). 

These trends show that the two-parameter TF relation does not adequately describe the 
data. We improve it by writing the TF relation in the following, four-parameter form: 

v = - e (M — D) — a(fx e - 19.2) + /3(c - 2) . (12) 

(We normalize fi e and c to typical values so as to minimize changes in the TF zero point.) 
The second and third rows of Table ^| show the results of incorporating first /i e and then 
both /i e and c into the TF relation. With each added free parameter, the likelihood is 
significantly improved (see the discussion in §5.1). In particular, incorporating /i e results 
in a 21 unit decrease of C, a J (21 — 1) ~ 4.5-sigma likelihood increase. Thus, the surface 
brightness dependence of the TF relation is highly statistically significant, although the 
maximum likelihood value of a ~ 0.05 is not especially large, corresponding to ^tf oc I®' 13 - 
The likelihood improvement when concentration index is incorporated into the fit is smaller, 
but still significant at the 3 a level. Note that this improvement occurs after SB has been 
added to the fit; thus, the SB- and concentration index- dependences of the TF relation are 
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Fig. 3. — Inverse TF residuals from a standard 2-parameter (rotation velocity-luminosity) TF fit 
are plotted against surface brightness fj, e (upper panel) and concentration index c (lower panel). 
The diagonal lines indicate the approximate slopes, 0.05 for fi e and 0.04 for c. The trends are 
significant at the ~ 4.5 a level for fi e and at the ~ 3a level for c (see text), indicating that the TF 
relation depends on both of these parameters. 

distinct from one another. Residuals with respect to /i e and c from the four-parameter TF 
relation, with the parameters given in the third row of Table |2|, are shown in Figure [|. 
No trends can be seen, indicating that equation [T2| is an acceptable description of the TF 
relation. Note that although f s was treated as a free parameter when fitting equation [12], 
its maximum likelihood value was virtually unchanged from the two-parameter fit (compare 
the values of f s listed in rows 1-3 of Table |2|). 
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Fig. 4. — Same as the previous figure, except that now the residuals for a four-parameter TF fit 



(equation 12) are shown. No trends are evident. 

5.7. Luminosity-Dependence of the TF Scatter 

Fitting the four-parameter TF relation involved six free parameters: the four in Eq. |i~2|, 
plus <j v and f s . A significant improvement in the likelihood can be achieved with the addition 
of a seventh parameter that describes the decrease in TF scatter for higher luminosity, higher 
surface brightness galaxies. Such a scatter decrease was detected for certain low-redshift TF 
samples by Giovanelli (1996), Willick et al. (1997a,b) and Willick & Strauss (1998). Here, 
however, the effect is even stronger, as can be seen in Figure |5|, in which TF residuals are 
plotted versus luminosity for the six-parameter fit. Low- luminosity galaxies exhibit larger 
residuals, in the mean, than do high-luminosity objects. We follow Willick & Strauss (1998) 
and write a v = cx^o + 9i(M — M) where M is the mean absolute magnitude of the sample 
(—21.48 for the LP10K sample). When g>j was included as the seventh free parameter in 
the fit, a 13.2-unit decrease in £ relative to the six-parameter fit was achieved, a 3.2a 
improvement. The best-fit result was = 0.011, more than twice as large as the value found 
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Fig. 5. — Inverse TF residuals from the four-parameter TF fit, equation 12, plotted versus lumi- 
nosity. It can be seen that low-luminosity galaxies exhibit larger TF scatter than high-luminosity 
ones. 



by Strauss & Willick (1998) for the MAT sample. Thus, TF scatter decreases by ~20% for 
each magnitude of luminosity. 

However, luminosity alone may not be the sole indicator of TF scatter. For example, in 
the upper panel of Figure | we see that scatter also correlates with surface brightness, in the 
sense that higher-SB galaxies appear to have smaller TF scatter. We can describe this effect 
with a linear model as well, and find a similar likelihood increase as we did in the previous 
paragraph. 

TF scatter thus decreases with increasing luminosity and with increasing surface bright- 
ness. Each effect alone is about equally significant. This suggests that the underlying factor 
is that the scatter decrease is most strongly correlated with increasing vtf itself. Such an 
effect indeed has a natural physical explanation. Suppose that one makes, on average, a 
constant error 5v^p in measuring the rotation velocity (including errors in the inclination 
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TABLE 2 

TF FIT PARAMETERS: / s -FORMULATION a,b 



D 


e 


a 


P 


/. 


a v 




C d 


Notes 


-21.619 


0.1314 






1.994 


0.0659 




-1068.1 


e,f,g 


-21.657 


0.1216 


0.0467 




2.033 


0.0631 




-1089.3 


f,g 


-21.765 


0.1102 


0.0533 


0.0454 


1.957 


0.0617 




-1100.0 


g 


-21.758 


0.1122 


0.0589 


0.0450 


2.001 


0.0366 


15.95 


-1116.8 


h 



Table 2: Notes: (a) The parameters D, e, a, and (3 are defined by equation [T^. a v is the overall 
inverse TF scatter, except as noted in the fourth row of the table, (b) The TF rotation velocity is 
defined according to equation |j| i.e., it is the rotation velocity evaluated from the arctan fit to the 
RC at / s times the effective exponential scale radius r e . f s is treated as a free parameter, and its 
maximum-likelihood value is listed, (c) The mean rotation velocity measurement error, in kms , 
as determined from maximum likelihood via equation |l4|. (d) The likelihood statistic is scaled by 
the ratio of the number of independent (245) to total (352) data points in the fit (see §5.1). (e) 
SB dependence of the TF relation not modeled, (f) Concentration-index dependence of the TF 
relation not modeled, (g) treated as a constant in fit. (h) Overall TF scatter modeled according 



to equation 14. The listed value of a v actually represents the quantity ctq in that equation. 



correction). The corresponding error in the circular velocity parameter is 

1 5v TF „„ rriri _„ ^TF 



5rj 



0.05510" 



20 km s- 



(13) 



In 10 v tf 

where in the last step the definition of rj has been explicitly taken into account. A reasonable 
model is then to assume that the total TF scatter consists of an ^-independent portion 
(intrinsic scatter plus photometric and distance assignment errors), o"o plus the error specified 
by equation O, added in quadrature: 



0V 



0n + 



0.055 



Sv' 



TF 



20 



(14) 



(Note that in evaluating a v one uses not the observed but the predicted r] in the exponent, 
as we are doing an inverse fit.) 

A fit using equation |14] for a n produced a greater likelihood increase than either a 
luminosity- or an SB-dependent scatter alone. The maximum-likelihood parameters resulting 
from such a fit are given in the fourth row of Table ||. The likelihood statistic changed by 
about 17 units with the addition of a single free parameter, 5v^f, a four-sigma improvement. 
Combined with its having a reasonable physical motivation, this likelihood increase leads us 



to adopt equation 14 as a meaningful description of the TF scatter. 
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Fig. 6. — Inverse TF residuals plotted against the predicted TF velocity, t>TF, based on the four 
parameter fit. The solid lines indicate the la, and the dashed lines the 3 a, predictions based on 
equation 14. With the exception of a larger-than expected number of 3 a deviations, equation 14 
gives a good account of the observed TF scatter. See text for details. 



The maximum likelihood values of the scatter parameters are <jq = 0.036 (corresponding 
to ~ 0.30 mag), and 8vtf = 15.6 km s _1 . The latter is a reasonable value for a typical rotation 
velocity measurement error. Figure ^| shows TF residuals plotted versus v^f (predicted) = 
158.1 x lov(M,n e ,c) kms -i_ Overplotted are the values of a v and 3 predicted from equa- 
tion [14|. The model provides a reasonable description of the observed TF residuals. The 
one way in which it breaks down is that there are more 3— a deviant points than one would 
expect from Gaussian statistics. This may point less to a failure of this particular model of 
the scatter, but instead to the general non-Gaussianity seen in many TF samples (Willick 
et al. 1996), in the sense that a few percent of any given sample is likely to exhibit large 
residuals. The remainder of the points are, to a good approximation, Gaussian. 

The curves in Figure ^ show that velocity width measurement errors dominate the 
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TF scatter for slow rotators (v^f ~ 150 kms -1 ), while for rapid rotators, i>tf ~ 200 kms -1 , 
i.e., high- luminosity, high SB galaxies, the constant error term, ero, dominates. This error, 
as noted above, includes photometric and distance errors, plus intrinsic scatter. Because 
the former are quite small in general (^0.1 mag), the value of o"o is largely indicative of 
the intrinsic scatter of the TF relation. Our result here thus implies that the intrinsic TF 
scatter is ;$ 0.30 mag, a result obtained through entirely independent means by Willick et al. 
(1996). The overall scatter within the observed range of vtf ranges from about <j v ~ 0.085 
for -utf ~ 100 kms -1 to <j v = <7o = 0.037 for utf ~300 kms -1 . Given the inverse TF slope 
of e = 0.112, these values correspond to equivalent forward TF scatter values of 0.33-0.75 
mag. This large range helps explain why different samples have yielded different values of 
the TF scatter in the past (cf. the discussion by Willick 1998c). 



5.8. An Alternative Formulation of the TF Relation 



The TF relation, equation [12], explicitly contains four parameters (e, D, a, /?), but 
implicitly a fifth parameter, f s , is present in the definition of wtf and thus rj. An equivalent 
approach is to make the fifth parameter explicit, as follows. First, we let r] a = log(2t> a ) — 2.5, 
i.e., we take the asymptotic velocity from the arctan fit as the TF rotation velocity. Then 
we make x t an explicit parameter in the TF relation: 

Va = -e(M - D) - a(fi e - 19.2) + (3{c - 2) + 7 x t . (15) 

To fit this TF relation to the LP10K data set we must first prune the sample a bit further 
by requiring x t < 2. For x t > 2 the RCs are very nearly linear and v a is highly uncertain 
(recall that purely linear RCs had already been eliminated). The resultant sample consists 
of 341 data points comprising 237 distinct galaxies. 

Figure [7| shows TF residuals [now defined as ^(observed) — r\ a (predicted)] from equa- 



tion y as a function of x t . In the top panel, residuals resulting from a four-parameter fit, 
with 7 = 0, are shown, while in the bottom panel residuals from the five-parameter fit (with 
the best-fit value of 7 indicated) are plotted. In each case, the scatter model given by equa- 



tion IT3 was used, with 5wtf held fixed at 15.95 kms 1 to facilitate comparison with the 



/ s -formulation. The maximum-likelihood parameters for the two fits are listed in Table |3|. A 
strong trend with x t is obvious in the upper plot; the lower plot shows no trend, indicating 
that a the linear model is a good one. The statistic £ decreases by 132 units when x t is 
included in the fit, a huge likelihood increase. 

These results show that when x t is not taken into account (7 = 0), luminosity is a 
poor predictor of v a . Galaxies with slowly rising rotation curves (x t ^ 1) have a much larger 
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Fig. 7. — Inverse TF residuals from the five-parameter TF fit, equation [l5|, plotted as a function 
of the RC shape parameter xt- In the top panel the coefficient of xt in the fit, 7, has been fixed 
at zero. When this is done, a strong trend of residuals with xt is evident, and the corresponding 
TF scatter is very large (see Table 3). In the bottom panel, 7 was varied to produce the maximum 
likelihood fit. No trend of the residuals is now evident, indicating that when xt is included in the 
TF fit, the asymptotic RC amplitude v a can be effectively used in the TF relation. 

asymptotic RC amplitude than their luminous mass suggests. The term 0.165xt corrects this 
deficit: luminosity (along with SB and concentration), in combination with x t , predicts the 
asymptotic RC amplitude well. Indeed, the TF scatter and likelihood for the fit based on 
Eq. [15| are nominally better than (though statistically equivalent to) those obtained from a 
fit based on a Eq. [12], the /^-formulation, for the same 341-galaxy sample. 

The explicitly five-parameter TF relation is useful because it clarifies why a TF relation 
based on v a and luminosity (with or without SB and concentration) is a very poor one: galax- 
ies of a given luminosity can have very different RC shapes, as determined by Xt- Figure [S], 
in which Xt is plotted against absolute magnitude (left panel) and surface brightness (right 
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TABLE 3 

TF FIT PARAMETERS: ar r FORMULATION a ' 6 



D 




e 


a 


P 


7 




<5t>TF 


C d 


Notes 


-21.E 


*34 


0.0881 


0.0767 


0.0104 




0.0650 


15.95 


-958.3 


e,f 


-21i 


m 


0.1173 


0.0551 


0.0406 


0.165 


0.0347 


15.95 


-1090.7 


f 



Table 3: Notes: (a) The parameters D, e, a, f3, and 7 are defined by equation [l5|. (b) The TF 
rotation velocity is taken to be the asymptotic value, v a . The RC shape parameter x% now enters 
explicitly into the TF relation, equation [l5|. (c) The contribution to from instrinsic scatter 
and photometric errors. The remainder comes from the 5vtf term according to equation |TJ|. (d) 
The likelihood statistic is scaled by the ratio of the number of independent (237) to total (341) 
data points in the fit (see §5.1). (e) The a^-dependence Q f the TF relation is not modeled, i.e., 7 
is held fixed at zero, (f) To facilitate comparison with the /^-formulation, <5utf is held fixed at 
15.95 kms -1 . 

panel), further underscores this point. Although there is a correlation of x t with luminosity, 
it is not a tight one. There is, perhaps surprisingly, no correlation between RC shape and 
surface brightness at all. Galaxies with given luminous properties can have a wide range 
of RC shapes, and vice versa. Unless the RC shape is accounted for, via the parameter x t , 
luminosity cannot be expected to predict the RC amplitude at large radii. 



5.9. A Radius-based TF Relation: the Spiral FP 

The multiparameter TF relations discussed above are reminiscent of the FP for elliptical 
galaxies (cf. §2). The latter, however, usually take the defining structural parameters to be 
galaxy size, surface brightness, and velocity width; here we have stayed closer to the TF 
relation and used luminosity in preference to size. These formulations are not independent, 
since luminosity oc l e r 2 e . Thus, our TF relation i>tf ^ L°' 28 I® AA implies t>TF oc I® A2 r° e m . 
Solving for r e this yields r e oc v^p 9 I~°' 75 . One may compare this with the elliptical FP 
relation (e.g., Jorgensen et al. 1996) r e oc a -1-24 /" 0-77 . The velocity width dependence is 
significantly stronger for the spirals, but the surface brightness dependence very similar. 

The above result assumes L oc I e T 2 , which may not be precisely true. In contrast to 
the procedure generally used with ellipticals, in the LP10K data set luminosity is derived 
independently from I e and r e . Furthermore, the multiparameter TF fits above included the 
variables f s and c, and it is not a priori obvious that they should enter the fTF-fV-fe relation 
in the same way. It is thus worthwhile to fit a standard FP relation directly to the data. We 
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Fig. 8. — The RC shape parameter xt plotted against absolute magnitude (left) and effective surface 
brightness (right). There is a trend of Xt with luminosity, in the sense that higher-luminosity objects 
have smaller values of xt in the mean, but the trend is not at all a tight one. No meaningful trend 
of xt with surface brightness can be seen. 

do so by writing 

V = s (log Re—T) — afi e + (3c , (16) 



where rj is derived from vtf as in §5.3. In equation [16], R e is the absolute diameter, in kpc, 
corresponding to the observed angular diameter r e . Distance is derived from a pure Hubble 
flow model as before. 

The results of carrying out such a fit are presented in Table |j; for simplicity, the TF 
scatter was treated as a constant for the fit. The fit parameters imply a spiral FP of the form 
t>TF -Re' 51 -^e' 39 ) or R-e ^ v^p 6 !^ ' 76 . Again, the surface brightness exponent is quite similar 
to the elliptical case, but the velocity dependence is significantly stronger. It is noteworthy 
that the best fit value of f s is 2.05, identical within errors to what was obtained from the 
luminosity-based fit. We may conclude that the radius-based fit and the luminosity-based 
fit reflect roughly the same physical properties of the galaxies. However, the likelihood 
statistic is considerably worse, and the TF scatter higher, for the radius-based fit than 
for the luminosity-based fit (the appropriate comparison is with the third row of Table |2|). 
Evidently R e is not as good an indicator of rotation velocity as luminosity. This reduces its 
utility as a distance indicator. In future papers in this series, the luminosity fit will be used 
in preference to the radius-based fit. 
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TABLE 4 

TF FIT PARAMETERS: RADIUS-BASED FP a 

T s a | g g £ d 
0.537 0.513 0.157 0.076 2.055 0.0644 -1079.5 



Table 4: Notes: (a) The parameters T, s, a, and (3 are defined by equation 16, where R e is measured 
in kpc. (b) The parameter f s has the same meaning here as described in the notes to Table 2. (c) 
TF scatter treated as a constant, (d) The likelihood statistic is scaled by the ratio of the number 
of independent (245) to total (352) data points in the fit (see §5.1). 



6. Discussion 



In this paper several key features of the TF relation have been identified. The first 
is that the TF relation is optimized — its scatter is minimized — when wtf is taken as the 
amplitude of the rotation curve at about two disk scale lengths (cf. Figure 0). This would 
not have been the case if spiral RCs possessed "universal" properties, as has been advocated 
by Persic, Salucci, & Stel (1996; PSS96). Under the PSS96 hypothesis, it is not merely the 
amplitude but the entire RC that is specified by luminosity: v(r) = vq(L) f (x; L) , where 
x = r/R opt and R opt is a characteristic optical radius. Were this so, the radius at which one 
chooses to evaluate the RC amplitude should be irrelevant, provided it is the same multiple of 
the optical radiusQ for all galaxies. Here we have affirmed the conclusions of Courteau (1997), 
Kasen (1997), and Simon (1998), that there is something special about the RC amplitude 
at two disk scale lengths as regards the TF relation. Luminosity is a poor predictor of the 
value of the RC at very small and very large radii. 

This result implies that RCs are not universal. Galaxies of a given luminosity can have 
RCs of very different shapes. This is exemplified by Figure which shows that the RC 
shape parameter x± is only loosely correlated with luminosity (and not at all with surface 
brightness). High- luminosity galaxies are more likely to have classical, flat RCs (xt <C 1) 
than are low-luminosity ones, but some high-luminosity galaxies have slowly rising RCs. Low 
luminosity-spirals are more likely to have quasi-linear RCs (x t ^l), but some in fact have 
flat RCs. This is why one must evaluate TF rotation velocity, vtf, at a particular radius, 
defined by the surface brightness profile, to obtain a strong correlation with luminosity. 

A corollary is that the asymptotic or flat value of the RC, v a , is not well correlated with 



5 Whilc PSS96 use a radius R opt defined by a fixed fraction of enclosed light, we have used a characteristic 
radius r e here. However, the two are linearly related with relatively little scatter and thus the distinction is 
unimportant. 
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luminosity, as shown especially well by Figure 0. If the RC rises slowly within the luminous 
disk (x t £jl), the luminosity underpredicts v a . In order to use v a as the TF velocity, one 
must explicitly incorporate the shape parameter xt into the TF relation, as in Eq. [Tj| The 
resultant TF relation exhibits scatter as small as that obtained by taking -utf as the RC 
amplitude at two scale lengths. 

The fact that TF scatter is minimized when v^f is taken as the RC amplitude at two 
disk scale lengths has important implications for the dark halos of spiral galaxies. Let us 
assume that spirals consist of a thin disk of mass surface density S(r) = T,Qe~ r ^ Rd , and a 
spherical dark halo with a density profile p(r) = phx( r /Rh)- Here, ph is a generic "central 
density" of the halo (not necessarily the density at r = 0, which is infinite in some models), 
and Rh a generic scale length (we consider a specific model for the dimensionless density 
profile x{ r /Rh) below). One can then show (Dorris & Willick 1998) that the rotation velocity 
at radius r may be written 



VAX) 



AirGT.oRd k(x) + 7a 1 e(x/a) 



where: 



x = 



Ra 



X 



«0) = -j 



where / and K are Bessel functions of the first and second kinds, 

PoRl 



7 



a = 



Rh 

Rd 



and the function e(x/a) is defined by 

1 rv 

e {y) = — I x x( x ) dx . 

V JO 



(17) 

(18) 
(19) 

(20) 
(21) 

(22) 



Note that 7 is, to within factors of order unity, the ratio of halo to disk mass, while a 
measures the relative sizes of the two components. 



It is evident from Equation \L7\ that, if a and 7 had "universal" values — if dark halos were 
always the same mass and size relative to the luminous disks which inhabit them — circular 
velocity would be a function of x = r/R d alone. If that were the case one could take wtf to 
be the amplitude of the RC at any multiple of the disk scale length, without changing the 
accuracy of the TF relation. A plausible explanation for the TF scatter being minimized at 
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x ~ 2 is that the halo parameters po and Rh vary for disks of a given mass and scale length 
and thus that a and 7 are not constant. 

We can quantify this idea further as follows. Suppose that halos vary, but in a predicable 
way, such that 7 is a function of a (i.e., central density is a function of physical size, or mass), 
as expected from theories of hierarchical structure formation (e.g., Navarro, Frenk, & White 
1996, hereafter NFW). Then from Eq. [17|one can show that variation in lnv c due to changes 
5a in halo size is 



5 InvJx) 



2vl 



x 



din 7 
din a 



1 + 



dine 



dlnz 



z=x/a j 



5a 



a 



(23) 



where v\ = 4:7iGpoRf l e(x/a) is the square of the rotation speed due to the halo alone. If 
indeed the TF scatter minimizes at x = 2 because halo mass and size variations have no 
effect there, by Equation E3] it follows that 



din 7 
din a 



1 + 



dine 



dlnz 



(24) 



--2 1 a 



This represents a constraint on the properties of the dark halo. If this constraint is satisfied, 
halo mass and size variations will not affect the scatter of the TF relation when v tf is taken 
at r = 2Rd- On the other hand, if Equation |24| is satisfied, then at some other radius x the 



RC amplitude variation due to halo mass and size variations will, according to Equation g3 
be given by 
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As a concrete example, let us suppose that the halo mass distribution is given by the 
model proposed by NFW, 
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Rh 
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PnfwIH = Po 

Using pNFw( r ) m Equation ^2] we find 

dlne(z) z 2 /(l + z 2 ) 
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dlnz 



ln(l + z) - z/(l + z) 
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(27) 



Before using this formula, we must adopt a suitable value of a. The requirement that 
RCs be well-approximated by an arctangent shape constrains a to lie in the range ~ 7- 
12 for luminous disks within NFW halos (Dorris & Willick 1998). We adopt a = 10 
for this discussion; the conclusions are not sensitive to its precise value. We then find 
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that (dlne/dhi z) z= 2/ a — 0.78. If we consider only luminous disks of fixed mass and size, 
d In 7/0? In a = dhiMh/dhxRh. If we further assume ph oc R^ , the constraint Eq. ^implies 
Ph oc Rh 1 ' 2 , or, equivalently, ph oc M^ ' 7 . Thus, the TF scatter minimum when t> TF is taken 
at two scale lengths, along with several additional assumptions, implies a particular scaling 
of halo central density with halo mass. The indicated scaling is consistent with the predic- 
tions of hierarchical structure formation scenarios, although a rather flat power spectrum 
(n — 1.5) on galaxy halo scales would be required (see Figure 9 of NFW). 



If we then use Equation |27| in Equation [25|, take x = 5 (where we assume the RC is 



halo-dominated), and assume that typical variations in a are 50%, we find that 5lnv c (x = 
5) ~ 0.08, corresponding to Srj ~ 0.034. Such an error, added in quadrature with a minimum 
error a v ~ 0.06, yields a v = 0.069. This is quite comparable to what is shown in Figure ^ as 
we go from f s — 2 to / s = 5. Thus, reasonable variations in halo mass at fixed disk mass can 
account for the increase in TF scatter when i>tf is taken as the asymptotic RC amplitude. 

The preceding argument is illustrative only. The NFW profile may not describe low- 
surface brightness spirals well (Kravtsov et al. 1998), and we have neglected the effect of a 
bulge component altogether, which is certainly important at small x. Nonetheless, the main 
point of the argument remains valid: if spirals of a given mass and size can be found within 
dark halos of variable mass and size, this will, in general, contribute to the TF scatter when 
v tf is taken at an arbitrary radius. However, if the halos exhibit a well-defined mass-size 
relationship, as embodied in Equation |24], this contribution to the TF scatter will vanish 
when v tf is taken as the RC amplitude at two disk scale lengths. A detailed consideration 
of these issues will be presented in a forthcoming paper (Dorris & Willick 1998). 

A second important result of this paper has been that the TF relation is not a two 
parameter relation. The v^f-L relation exhibits residuals that correlate, though weakly, 
with surface brightness (upper panel of Figure [3]). When SB is incorporated into the relation 
one obtains a significant reduction in scatter and improvement in likelihood of the fit. The 
best-fit TF relation may be written v^f oc L 0,28 /^' 14 , where I e is effective central surface 
brightness (cf. Appendix A). If we neglect mass-to-light variations, (i.e., if we assume mass 
is proportional to luminosity and surface density is proportional to surface brightness) then 
Eq. |17] reduces to i>tf °c L°' 25 I® 25 . The observed TF relation is very close to this virial 
prediction in terms of the luminosity-dependence, but the surface brightness dependence 
is significantly weaker. Nonetheless, the existence of an SB-dependence is important. It 
has generally been assumed that the TF relation lacks such a dependence, and this has 
been explained in the past as a consequence of the dark halo dominating the mass within 2 
scale lengths and masking the SB dependence (Courteau & Rix 1999; McGaugh & de Blok 
1998, Navarro 1998) The detection here of a significant SB dependence of the TF relation 
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shows that the virial theorem is approximately, though not precisely, obeyed within two 
scale lengths. The halo may conspire to reduce the SB dependence, but it does not eliminate 
it altogether, which requires less fine-tuning to explain. Alternatively, the sub-virial SB 
dependence could be a consequence of a mild decrease in mass-to-light ratio with increasing 
surface brightness. Which if either of these possibilities holds is not clear; the problem 
deserves further study. 

A related conclusion we may draw from this result is that the properties of spirals are 
closer to those of ellipticals, which exhibit fundamental plane relations, than has previously 
been appreciated. This conclusion has been anticipated in the work of Bender, Burstein, 
Faber, & Nolthenius (1997), who have emphasized that all galaxies — indeed all gravi- 
tationally bound stellar systems — inhabit a fundamental plane, the "cosmic metaplane." 
However, Bender et al. considered the two-dimensional, vtf-L, locus in parameter space, 
to represent the FP of spirals. The SB-dependence of the TF relation shows that the FP 
of spirals is not parallel to the vtf-L plane. We have also shown that the spiral data may 
be fitted directly to the elliptical FP variables (radius, velocity, and surface brightness), 
yielding the result R e oc t^p 6 /" ' 76 . The I e exponent is very similar to what has been found 
for ellipticals, though the v^f exponent is significantly larger. However, we found that the 
FP defined by i>tf, logi? e , and I e (the "radius-based" FP) exhibits noticeably larger scatter 
than the v^F-L-I e fit (the "luminosity-based" FP). In this sense, spirals hew more closely 
to the traditional TF relation than to the FP of ellipticals. This may well be related to the 
fact that, because of their relatively sharp edges, it is easier to measure the total luminosity 
of spirals than it is for ellipticals. 

A final comment concerns the dependence on concentration index, c, found here for both 
the luminosity and radius-based spiral FPs. This is probably not a fundamental considera- 
tion. Rather, the presence of c probably serves mainly to correct the effective radius r e for 
galaxies of different luminosity profiles. Its effect is significant only for highly concentrated 
galaxies. These are objects with strong bulges for which r e , which is relatively insensitive to 
the innermost SB profile, overstates the spatial extent of luminous mass. The concentration 
index dependence of the TF relation represents a small correction for this effect. 



7. Summary 

We have presented the first results from the Las Campanas-Palomar 10,000 kms -1 
(LP10K) cluster survey. The survey includes TF data for spirals and FP data for ellipticals 
found within ~ 1 deg of the centers of 15 Abell clusters. The elliptical galaxy observations 
and analysis are ongoing; this paper, along with Papers II and III, concerns the TF results 
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only. The principal aim of the LP10K survey is to constrain bulk peculiar velocities on very 
large ( £ 100/i _1 Mpc) scales. 

The LP10K TF data set is one of a number of recent TF surveys based on rotation 
curves measured from long-slit spectroscopy. We may now ask, how can we use the full 
information contained in the RC — and not just a single "velocity width" — to optimize the 
TF relation? Another relevant question is how best to use detailed surface photometric 
information derived from CCD imaging in TF analyses. The focus of this paper has been to 
address these and related issues using the LP10K TF data set, and to interpret the results 
in light of plausible models of spiral galaxy structure. 

Several key parameters were derived from the photometric and spectroscopic data. In 
addition to an i?-band apparent magnitude, the surface brightness profile was used to cal- 
culate an effective exponential radius r e and central surface brightness I e . These parameters 
were determined from intensity moments (Appendix A) rather than from a direct fit to the 
profile, and thus are robust and objective. A luminosity concentration index c was also com- 
puted from the surface photometry profile. The long-slit spectroscopy produced a rotation 
curve v c (r) measured from the Ha emission line. The RCs were fitted by a two-parameter 
arctangent function (Equation [|), an adequate fit in all cases. This procedure yielded an 
RC turnover radius r t and asymptotic velocity v a . We define Xt = r t /r e , the ratio of the 
dynamical to luminous scale length of the galaxy. Objects with j; t < 1 exhibit classical 
rotation curves with a flat or nearly flat portion. Objects with x t ^ 1 have slowly rising RCs 
that may not flatten out within the measured region. 

We carried out a series of maximum-likelihood fits of the inverse TF relation, essentially 
minimizing residuals of the circular velocity parameter rj = logt>TF — 2.5. The TF rotation 
velocity t>TF was obtained from the arctangent fit: 



The quantity f s was treated as a free parameter to be determined via likelihood maximiza- 
tion. The fits were optimized for f s = 2.0 ± 0.2. This conclusion confirms that reached earlier 
by Courteau (1997), Kasen (1997), and Simon (1998). These studies indicate that vtf must 
be referenced to the luminosity structure of the galaxy, not taken from the asymptotically 
flat part of the RC; only when i t <Cl are v^f and v a effectively equivalent. This finding is 
at variance with the conventional wisdom concerning the TF relation, which holds that it 
is the amplitude of the flat portion of the RC which best correlates with luminosity (e.g., 
McGaugh & de Blok 1998). A corollary is that spirals cannot be characterized in terms of a 
"universal" rotation curve (URC), as has been argued by Persic, Salucci, & Stel (1996). The 
URC hypothesis holds that luminosity uniquely specifies not only RC amplitude, but also 
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the entire functional form of the RC. Were this the case, one should be able to evaluate utf 
and any chosen multiple of the effective radius r e , and obtain an equally tight TF relation. 
As Figure 2 clearly shows, this is not the case. 

The physical mechanism most likely responsible for the TF optimization at two scale 
lengths is variation in the mass and size of the dark halos in which galaxies of given luminous 
mass and size reside. A simple model was presented (§6) to demonstrate the outlines of this 
effect. If a spherical dark halo of central density p^ and scale length Rh surrounds a disk 
galaxy, variations in ph and Rh will affect its contribution to the RC at all radii. However, 
these effects will cancel at r = 2r e , provided that ph varies with Rh in the manner predicted 
by equation This constraint is dependent on the specific model of the halo one adopts. 
For an NFW halo whose parameters are tuned to yield reasonable RC shapes, the constraint 
may be stated pu oc M^ 0,7 , where Mh is the total mass of the halo. Such a result is roughly 
consistent with the mass-density relation derived from numerical simulations of structure 
formation (NFW). 

We have also demonstrated that TF residuals from a standard two-parameter fit, t>TF (X 
L a , are correlated, albeit weakly, with surface brightness. A tighter TF relation is obtained 
if SB is explicitly incorporated into the fit. The resultant three-parameter TF relation may 
be written vtf L 0,28 /g 14 . Such a TF relation resembles the fundamental plane of elliptical 
galaxies more closely than does the standard two-parameter relation, and indeed is closer to 
the prediction from simple virial equilibrium considerations. However, the surface-brightness 
dependence remains sub-virial, a result still in need of explanation. 

We found a significant decrease in the TF scatter with increasing luminosity and surface 
brightness. We suggested that this effect follows naturally if rotation velocity measurement 
errors, Sv-ff, are independent of v^f itself, so that <51nt> TF cx v^ F . Incorporating this model 
improves the fit likelihood and suggests that typical rotation velocity measurement errors 
are ~ 15 kms" 1 . These errors dominate the TF scatter for slow rotators (i>tf ~ 150 kms -1 ), 
while for rapid rotators — high luminosity, high-SB galaxies — the intrinsic scatter of the TF 
relation, ^0.30 mag, dominates the observed scatter. 

The above properties of the TF relation will be used in subsequent papers in this series, 
which will investigate deviations from uniform Hubble flow using the LP10K sample. 

The author acknowledges the support of NSF grant AST-9617188 and the Research 
Corporation, and thanks Michael Dorris, Stephane Courteau, Michael Strauss, Mike Hudson, 
and the anonymous referee for valuable comments on earlier drafts of this paper. 
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A. Characteristic Surface Brightness and Scale from Moments of the 

Intensity Distribution 



It is customary to model spiral galaxy intensity profiles as exponential, J(r) = I$e 



r/r d 



The central intensity Iq and the disk scale length then provide convenient measures of 
surface brightness and radius. For two reasons, however, it is not optimal to obtain these 
quantities by fitting exponential functions to the observed intensity profile. First, many 
real spirals are not accurately exponential. Second, any fitting of a profile is necessarily 
subjective: a start and end point for the fit must be selected, particular points may be 
rejected, etc. Different choices will yield different values of I and of r d . 

For these reasons, an objective measurement of characteristic surface brightness and 
scale — one that does not rely on a fit — is needed. However, given that many spirals are 
approximately exponential, it is desirable that the objective approach recover the exponential 
parameters in the event that the galaxy is, in fact, accurately exponential. When it is not, 
the approach should still produce well-defined and sensible measures of surface brightness 
and scale. Such a method is described in what follows. 

Given an intensity profile J(r) measured to a maximum radius r/, define the ra-th mo- 
ment of the intensity distribution by 



The moments f n are readily computed from the measured surface brightness profile, usually 
by taking I(r) as the median surface brightness measured on an elliptical isophote of major 
axis radius r. 

If I(r) is indeed exponential with central intensity Jo and scale length r<2, then the 
moments are given by 




(Al) 




(A2) 



where yf = rf/rj, and the dimensionless functions h n are defined by 




(A3) 



From equation |A2| we find that, for an exponential profile, 



}_h = h(y f ) = 1_ 6-e-yf(6 + 6y f + 3y 2 f + y 3 f ) 
rjfi hiyf) yj l- e -yf(l + y f ) 



(A4) 
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where in the last step the integrals have been explicitly evaluated. From equation S2 we 
also find that, for an exponential profile 

h 



(A5) 



Equation [A4] is a transcendental equation for yf in terms of the observables r/, fi, and 
/ 3 that may be solved using standard numerical techniques. Having obtained yf in this way, 
one calculates the exponential scale length = rj/yj. The central surface brightness Jo is 
then obtained by substituting yf into equation |A5| . 

This procedure has been tested on simulated galaxy profiles. When the profiles are 
exponential, it recovers accurately the true central intensity Iq and exponential scale length 
rd- However, the procedure is perfectly well-defined for well-behaved but otherwise arbitrary 
intensity profiles as well. In particular, one may always calculate yf from equation |A4| , 
and then define an effective radius r e = Tf/yf and an effective central surface brightness 
h = jVlT/^G//)]- These effective quantities are objective and robust against irregularities 
in the profile. Most importantly, they are, as discussed in the main body of the paper, found 
to be suitable measures of radius and surface brightness in the TF relation. 
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